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1  Introduction 


The  development  of  multi-axis  force  sensors,  i.e.  instruments  for  measuring  several  (up  to 
6)  components  of  force  and  torque  simultaneously,  has  been  initially  undertaken  in  fields 
such  as  wind-tunnel  testing,  adaptive  control  of  machines  and  thrust  stand  testing  of  rocket 
engines.  Typical  in  those  areas  is  the  problem  of  monitoring  forces  of  variable  directions 
and  intensity.  Some  of  the  most  interesting  sensors  designed  in  this  phase  arc  reviewed  by 
Doebeiin  [1],  Starting  from  the  mid  seventies,  a  major  impulse  to  research  on  such  sensors 
has  been  given  by  the  recognition  of  their  large  relevance  to  robotic  and  telemanipulation 
applications. 

The  ability  to  sense  the  arm-environment  interactions  is  a  crucial  need  for  robots  to  evolve 
from  purely  repetitive  behaviors  to  some  degree  of  autonomy  in  unstructured  surroundings 
or  tasks.  The  feed-back  of  forces  and  torques  exerted  by  the  arm’s  end  effector  is  not  only 
instrumental  for  the  accomplishment  of  most  tasks  involving  modification  of  the  environment 
by  contact,  but  is  also  critical  in  guaranteeing  safe  operation  of  the  arm. 

Considering  the  most  common  case  of  a  serial  link  manipulator,  a  complete  characteriza¬ 
tion  of  the  system  of  forces  acting  on  a  portion  of  the  arm  can  be  obtained  by  force  sensors 
interposed  between  that  portion  and  the  rest  of  the  arm.  Force  sensing  at  intermediate 
sections  of  the  kinematic  chain  is  likely  to  be  very  useful  for  next-generation  robots  using 
the  whole  arm  surface  to  interact  with  the  environment  [2,  3).  However,  since  the  measure¬ 
ment  of  interaction  forces  can  be  disturbed  by  inertial  forces  arising  from  acceleration  of  the 
masses  between  the  sensor  and  the  end-effector,  force  sensors  are  usually  placed  as  close  as 
possible  to  the  distal  end  of  the  arm,  most  often  at  the  wrist-  Several  configurations  have 
been  proposed  for  such  sensorized  wrists  (e.g.  [4,  5,  6]),  and  some  are  commercially  available. 
A  pedestal  force  sensor  has  also  been  described  in  [4]. 

The  increasingly  demanding  tasks  assigned  to  automatic  manipulation  require  ever  more 
stringent  accuracy,  stiffness,  encumbrance,  weight,  speed  of  response,  ease  of  interfacing  and 
reliability  of  force  sensors  for  robotics.  A  particularly  promising  application  of  force  sensing 
to  fine  manipulation,  for  instance,  is  the  so-called  intrinsic  contact  sensing  concept  [7,  8]. 
Very  detailed  information  about  the  contact  between  two  surfaces  can  be  obtained  based  on 
force  measurements  and  geometric  considerations.  By  integrating  force  sensors  in  the  very 
“fingertips”  of  the  arm,  interactions  with  manipulated  objects  can  be  monitored  very  closely. 
The  improvement  in  manipulation  dexterity  has  proven  considerable  for  both  simple  parallel 
jaw  grippers  [9],  and  for  dextrous  articulated  hands  [10,  11,  12].  A  drastic  reduction  of  size 
and  weight  of  the  force  sensors  with  respect  to  existing  designs  turned  out  to  be  mandatory 
for  this  application. 

The  stringent  and  conflicting  requirements  of  this  and  other  advanced  applications  of 
force  sensing  render  unsatisfactory  the  traditional,  intuitive  approach  to  force  sensor  design, 
and  motivated  the  investigation  of  a  more  systematic  approach.  In  this  paper  the  design 
of  multi-axis  force  sensors  is  considered  within  the  framework  of  optimal  design  theory, 
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the  branch  of  operations  research  currently  offering  the  most  developed  tool  for  systematic 
synthesis  of  optimized  designs. 


2  Problei  i  Formulation 

According  to  Vanderplaats  [13],  an  optimal  design  problem  (which  is  basically  a  nonlinear 
programming  problem)  can  be  stated  as  follows: 

Find  the  set  of  r  design  variables  x;  that  minimizes  the  objective  function  F(x), 
subject  to: 

1.  p  inequality  constraints:  <?i(x)  <  0,  1  <  i  <  p; 

2.  q  absolute  constraints:  hj(x)  =  o,  i  <  i  <  ?: 

3.  n  bounds  defining  the  feasible  design  region:  !*,,„/  <  x*  <  x*i<up,  1  <  k  <  n. 

The  variables  characterizing  a  force  sensor  design  can  be  chosen  within  a  large  range  of 
quantities  related  to  its  geometric  description,  dimensions,  material  properties,  machining 
process,  transduction  principle,  etc..  Any  quantity  contributing  to  the  identification  of  a 
design,  which  is  not  allowed  to  vary  in  the  optimization  process,  is  considered  a  design 
parameter. 

Inequality  conditions  are  unilateral  constraints  that  must  be  satisfied  in  order  for  the 
design  to  be  acceptable.  For  example,  stresses  in  structure  members  must  not  exceed  specified 
values.  The  maximum  displacement  under  nominal  load,  and  the  resolution  of  the  sensor 
(intended  here  as  the  the  minimum  detectable  load)  can  be  conveniently  considered  under 
the  form  of  inequality  constraints. 

By  means  of  an  equality  constraint,  precise  conditions  that  must  be  met  by  the  design 
can  be  stated.  For  instance,  these  constraints  can  be  used  to  define  the  system  interfaces.  A 
complex  problem  can  be  reduced  to  a  simpler  sub-optimal  one  by  using  equality  constraints 
to  assign  some  variables  a  fixed  value. 

Finally,  boundary  constraints  defining  the  feasible  region  for  design  variables  prevent 
the  optimization  algorithm  from  converging  to  unacceptable  solutions  (e.g.,  sensor  structure 
members  with  negative  thickness). 

The  objective  function  F(x)  is  the  measure  adopted  for  quantifying  the  quality  of  a 
design.  A  most  important  decision  is  implied  by  the  choice  of  the  quality  to  be  optimized. 
In  some  cases  emphasis  can  be  put  on  cost,  encumbrance,  or  weight  reduction.  Very  often, 
however,  the  accuracy  of  measurements  is  of  paramount  importance  in  designing  a  force 
sensor.  It  is  possible  to  take  into  account  these  and  other  possible  criteria  at  the  same  time 
by  using  multi-criteria  optimization  techniques  [14],  e.g.  by  assigning  weighting  factors  to 
each  criterion.  Although  greater  flexibility  of  the  optimization  tool  can  be  achieved  in  this 
way,  there  is  a  risk  with  multi-criteria  techniques  of  loosing  insight  into  the  design  process.  In 
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this  paper,  ail  the  above  mentioned  and  other  possible  requirements  are  considered  as  design 
constraints,  and  the  focus  is  on  the  choice  of  an  objective  function  capable  of  extracting, 
from  the  class  of  acceptable  designs,  the  one  giving  best  results  in  terms  of  measurement 
accuracy.  In  order  to  obtain  general  solutions  from  the  optimal  design  method,  we  propose 
an  objective  function  that  covers  a  wide  scope  of  design  variables.  In  particular,  an  attempt 
is  made  at  addressing  the  choice  of  the  number  of  basic  transducers  to  be  employed  in  a 
multi-component  sensor. 

3  Mathematical  Model  of  a  Force  Sensor 

In  general,  a  multi-axis  force  sensor  is  a  device  in  which  several  simple  transducers  measure 
the  effects  of  unknown  loads  on  a  mechanical  structure.  These  basic  measurements  are 
processed  in  order  to  evaluate  the  components  of  the  applied  load.  As  far  as  a  linear  behavior 
can  be  hypothesized  for  the  sensor,  a  model  can  be  written  as 


V  =  Cp  (1) 

where  p  €  3s?” ,  n  <  6,  is  the  unknown  load  vector,  V  6  3?m,  m  >  n  is  a  vector  collecting 
the  basic  measurements,  and  C  G  f?mXn  is  a  constant,  full  rank  matrix  characteristic  of  the 
sensor.  The  unknown  vector  p  is  composed  of  some  or  all  the  components  of  the  resultant  of 
the  load;  each  component  is  normalized  with  respect  to  its  maximum  value,  which  is  given 
as  a  design  specification.  When  basic  measurements  are  relaiive  to  strains  or  displacements 
(as  it  most  often  occurs  in  practical  applications),  the  C  matrix  is  usually  referred  to  as  the 
compliance  matrix  of  the  sensor. 

The  behavior  assumed  in  equation  1  is  ideal  under  several  regards.  For  instance,  the 
linearity  assumption  is  rather  strong.  However,  it  should  be  pointed  out  that  we  are  forced 
to  such  approximation  by  the  need  for  a  viable  algorithm  for  inverting  the  model,  i.e.  for 
solving  equatioh  1.  Nonlinear  elastic  models  of  even  simple  mechanical  structures  are  too 
complex  for  being  inverted,  and  even  though  a  tabulated  calibration  of  a  sensor  is  conceivable, 
there  is  no  practical  application  of  such  technique  to  force  sensors.  As  a  matter  of  fact,  linear 
(or  piece-wise  linear)  approximations  can  be  sharpened  by  calibration  techniques  (see  e.g. 
[15])- 

As  a  second  remark  about  the  model  1,  we  note  that  the  measurement  vector  v  will  be 
in  all  likelihood  affected  by  an  error  6v  due  to  inaccuracies  of  basic  transducers.  Similarly, 
also  the  imprecise  knowledge  of  the  actual  C  matrix,  which  will  be  affected  by  a  (matrix) 
error  SC,  should  be  taken  into  account.  Thus,  a  model  for  a  real  sensor  should  be  written  as 


v  +  8v  =  (C  +  SC)(p  +  Sp),  (2) 

where  £p  is  the  error  resulting  on  the  ultimate  information  from  the  sensor,  the  load  vector 
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Consider  for  instance  the  origin  of  errors  6v  and  SC  for  a  strain  gauge  based  sensor. 
As  is  well  known,  strain  gauges  are  electrical  resistors  that  can  be  bonded  on  a  deformable 
structure,  and  vary  their  resistance  according  to  the  mechanical  strain  of  the  structure. 
Errors  Sv  in  evaluation  of  structure  strains  can  be  caused  by  electrical  noise  and  thermal 
drift  in  resistance  measurements,  and  by  discretization  errors  in  conversion  of  data  from 
analog  to  digital.  The  elements  of  the  compliance  matrix  C  are  in  turn  determined  only 
approximately,  either  using  elasticity  theory  formulas  or  a  direct  calibration  procedure.  Both 
these  methods  result  in  a  relative  error  ec  (usually  much  larger  in  the  first  case),  arising 
from  inappropriate  modeling  of  sensor  structure  (for  calculated  C’s),  or  from  inaccuracies 
in  calibration.  The  effects  of  nonlinearities  in  the  stress- strain  relationship  of  the  structure, 
the  imperfect  stiffness  of  the  glue  used  to  bond  the  gauges,  and  many  ocher  uncontrollable 
factors  will  also  add  to  ec. 

Multi-axis  force  sensors  can  be  subdivided  between  those  using  the  minimal  number  of 
basic  transducers  for  the  measure  of  the  unknown  load  components  (i.e.,  m  —  n),  and  those 
having  redundant  sensors  (m  >  n).  In  the  latter  case,  no  exact  solution  of  equation  1  is 
possible,  and  diverse  approximate  solutions  can  be  chosen.  This  classification  will  be  used 
in  the  following  discussion. 

4  Minimal  Sensors 

For  rn  =  n,  the  generalized  form  of  Wilkinson’s  formula  for  error  propagation  in  linear 
algebraic  systems  (see  [16])  can  be  utilized  in  order  to  give  an  a  priori  estimate  of  the 
relative  error  on  p: 


ep  =  (e„  +  ec)  Kp( C),  (3) 

where  ^  ec  =  ep  =  are  respectively  the  relative  errors  on  strain  measure¬ 

ments,  on  calibration  and  on  the  results.  The  symbol  [j  -  ||  represents  a  generic  vector  norm, 
or  the  subordinate  matrix  norm  if  the  argument  is  a  matrix  1.  The  propagation  factor  Kp 
has  the  upper  bound: 


*p(C)  < 


Nc 


1-iW 

where  Nc  is  the  condition  number  of  the  compliance  matrix  C,  defined  as 


(4) 


Nc  =  II C}|  lie-1!!, 


and  where  it  is  assumed  Nc  ec  <  1. 

‘Vectors  are  designed  by  lower  case  boldface  letters  in  this  paper,  while  matrices  are  upper  case  boldface 
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A  block  scheme  of  equation  3  is  given  in  figure  1.  It  appears  from  equation  4  that  large 
condition  numbers  of  C  can  spoil  even  the  most  accurately  measured  and  calibrated  sensor. 

After  its  definition,  Nc  >  1,  which  implies  Kp  >  1:  this  can  be  seen  as  an  instance  of  the 
general  principle  of  information  theory  regarding  the  entropy  increr^se  due  to  elaboration  of 
data  [17].  In  the  best  case,  when  Nc  =  1  (which  occurs  only  for  orthogonal  matrices  and 
their  multiplies),  and  Kp  =  (1  —  ec)-1  ^  1,  the  global  error  ep  is  simply  the  sum  of  source 
errors  ev  and  ec. 

Given  a  sensor  design,  substantial  reductions  of  e„  and  ec  can  be  achieved  only  by  using 
more  sophisticated  technologies,  materials  and  components  in  the  construction  of  the  sensor, 
and  finer  models  of  the  structure,  or  more  accurate  instrumentation  for  strain  measurement 
and  calibration.  An  improvement  of  these  factors  with  respect  to  present  sensors  is  therefore 
achievable  by  increasing  their  cost.  Also  assuming  no  cost  constraint  to  the  design,  how¬ 
ever,  there  are  absolute  upper  bounds  to  possible  reduction  of  source  errors  set  by  present 
technological  state-of-art,  and  by  inherent  measurement  accuracy  limitations. 

The  amplification  factor  Kp  in  equation  3  is  related  to  the  elements  of  the  compliance 
matrix  C,  which  can  be  chosen  by  the  designer  by  varying,  for  instance,  the  thickness  of 
some  members  in  the  structure,  the  position  of  strain-gauges,  etc.  Thus,  the  problem  of 
optimizing  the  accuracy  of  multi-axis  force  sensors  can  be  split  in  a  technology-dominated 
sub-problem,  and  a  design-dominated  one. 

It  should  be  pointed  out  that  equation  3  does  not  assume  any  specific  algorithm  for 
the  solution  of  the  linear  system  1.  From  this  viewpoint,  different  algorithms  such  as  LU 
factorization,  singular  value  decomposition  (SVD),  and  even  special  solution  schemes  relying 
on  a  possibly  structured  C  matrix,  are  equivalent.  The  application  of  these  methods  with 
finite  precision  machines  will  produce  numerical  errors  that  are  typically  some  orders  of 
magnitude  smaller  than  source  errors,  and  can  be  considered  as  a  small  portion  of  the 
overall  calibration  error  ec  [18]. 

Since  the  right-hand  member  of  equation  4  is  a  monotonic  non-decreasing  function  of  Nc, 
it  seems  natural  to  choose  the  condition  number  of  the  compliance  matrix  C  as  the  objective 
function  of  the  optimization  procedure.  This  choice  is  indeed  correct  when  only  minimal 
sensors  are  taken  into  consideration.  M.  Uchiyama  and  K.  Hakomori  proposed  the  use  of 
the  condition  number  of  the  compliance  matrix  of  a  force  sensor  in  [19].  Unfortunately,  their 
work  was  published  in  Japanese,  and  it  has  not  been  as  widely  known  as  it  deserved  until 
cited  by  [20]  2.  The  authors  of  [20]  criticize  the  condition  number  criterion  of  [19]  under  four 
regards.  Some  of  those  points  are  discussed  here,  since  this  is  believed  to  provide  insight  in 
the  general  problem  of  optimal  design  of  force  sensors. 

The  first  remark  is  that  the  condition  number  criterion  is  insensitive  to  absolute  values 
of  the  compliance  matrix  entries.  In  other  words,  N(C)  =  N( 2C).  While  it  is  an  obvious 
advantage  of  condition  numbers  to  be  scale-independent,  the  authors  of  [20]  underscore  the 

’An  independently  developed  formulation  of  the  method  has  been  presented  at  the  same  time  in  [11] 
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fact  that  sensors  with  larger  absolute  values  of  structural  strains  are  not  rewarded  by  this 
criterion.  A  design  principle  consisting  in  making  the  “strain-gauge  sensitivity”  as  large  as 
possible  is  proposed,  which  consists  in  maximizing  the  norm  of  the  rows  of  C.  In  other  words, 
a  sensor  is  optimal  under  this  principle  if,  for  the  maximum  nominal  load,  every  gauge  is 
strained  at  the  maximum  allowable  level. 

The  second,  third  and  fourth  remarks  to  the  condition  number  criterion  all  apply  to  the 
fact  that  it  does  not  generalize  to  redundant  sensors.  For  instance,  the  two  compliance  ma¬ 


trices  Ci 


/  1  0 
1  0 

V  o  i  / 


are  brought  to  our  consideration.  The  condition 


numbers  are  Ac,  =  1  and  Nct  —  respectively.  The  authors  note  that,  “although  Cj  may 
have  a  redundant  strain-gauge,  C2  seems  not  to  have  less  sensing  performance  than  CV’. 
This  fact  is  true  indeed,  and  it  will  be  shown  in  the  following  sections  that  this  problem 
(as  well  as  the  others  pointed  out  by  [20])  can  be  solved  by  a  correct  generalization  of  the 
condition  number  criterion  to  redundant  sensors. 

Two  more  design  principles  are  proposed  by  the  authors  of  [20],  the  first  being  to  maximize 
the  minimum  singular  value  of  the  compliance  matrix  C,  and  the  second  to  minimize  the 
maximum  singular  value  of  a  displacement  matrix  G,  evaluated  at  the  point  of  application 
of  the  load. 

In  the  method  presented  in  this  paper  we  avoid  the  definition  of  multiple  optimality 
criteria,  in  order  to  avoid  the  need  of  a  difficult  choice  of  weights,  which  always  imply  a  degree 
of  arbitrariness.  A  lower  bound  on  sensitivity,  as  well  as  an  upper  bound  on  displacement 
under  load,  can  be  very  easily  incorporated  as  unilateral  constraints  to  the  design  (see 
section  2  above  and  7  below).  On  the  other  hand,  an  attempt  is  made  at  generalizing  the 
condition  number  criterion  for  redundant  sensors,  as  the  single  criterion  able  to  synthesize 
accuracy  requirements  for  the  sensor  (it  can  be  seen  for  instance  that  maximizing  strain 
gauge  and  force  sensitivities  as  defined  by  [20]  leads  to  minimizing  the  condition  number). 


5  Redundant  Sensors 

Consider  again  equation  1,  v  =  Cp,  in  the  case  that  the  number  of  independent  measure¬ 
ments,  m,  is  larger  than  the  number  of  unknown  load  components,  n.  Since  the  entries  of 
both  C  and  v  derive  from  experimental  measurements,  equation  1  is  in  any  practical  case 
inconsistent,  i.e.  no  exact  solution  is  possible.  Only  approximate  solutions,  denoted  with  p, 
can  be  obtained  through  a  manipulation  of  equation  1  as: 


=  MCp, 

=  (MCJ-'Mv, 


Mv 

P 


t 


(5) 

(6) 


Figure  1:  Error  propagation  in  a  force/torque  sensor 

where  M  6  $?nxm,  and  it  is  assumed  that  det(MC)  ^  0. 

Such  a  manipulation  consists  of  projecting  the  measured  point  v  of  the  m- dimensional 
space  onto  the  range  space  of  C,  denoted  with  7£(C),  which  is  a  n-dimensional  subspace.  The 
direction  along  which  the  projection  is  made  is  specified  by  the  choice  of  M.  In  particular, 
if  M  -  Cr,  the  projection  results  normal  to  7£(C),  and  the  corresponding  solution  p  = 
C^v  (where  =  (CrC)-1CT  is  th**  pseudo-inverse  of  C)  is  the  least  squares  solution 
to  equation  1.  Note  that  we  used  here  the  assumption  that  C  is  full  rank,  which  is  not 
restrictive  at  all  insofar  as  reasonable  sensor  designs  are  considered. 

Naturally,  any  minimal  sensor  can  be  considered  as  a  particular  case  of  a  redundant  one, 
by  supposing  that  (m  -  n)  measurements  are  neglected.  Correspondingly,  and  modulo  a 
reordering  of  rows,  a  manipulation  matrix  M  =  (InxnlOnx(m-n))  is  adopted  in  equation  5. 
Also  for  minimal  sensors,  therefore,  the  solution  p  is  to  be  considered  only  an  approximation 
of  the  “true”  value  of  applied  load. 

The  total  error  <5p  can  hence  be  thought  as  a  sum  of  two  terms,  deriving  from  the 
projective  manipulation  5,  and  from  the  actual  solution  of  6,  respectively.  The  rest  of 
this  paper  will  discuss  the  application  of  the  generalized  condition  number  criterion  to  the 
solution  of  equation  6:  minimizing  the  error  in  its  solution  is  assumed  as  the  goal  of  design 
optimization.  Before  beginning  such  discussion,  however,  it  seems  appropriate  to  show  that 
the  error  term  deriving  from  the  projection  5  is  not  directly  affected  by  the  choice  of  the 
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Figure  2:  The  projection  error  is  not  affected  by  the  design  of  the  compliance  matrix. 

compliance  matrix  C  elements,  but  rather  it  depends  only  on  the  source  errors,  <$v  and 
6C.  In  other  words,  the  projection  error  is  technology-dominated,  while  the  solution  error 
is  design-dominated. 

This  fact  can  be  graphically  illustrated  with  reference  to  a  very  simple  example  with 
m  =  2,  n  =  1.  An  hypothetical  error-free  sensor,  modeled  by  the  iinear  equation  v  =  Cp. 
would  have  the  measurement  vector  v  lying  exactly  on  'R(C),  represented  in  figure  2  with 
a  dashed  line.  The  “true”  solution  would  be  represented  in  this  case  by  the  length  of  the 
segment  OZ.  However,  any  solution  of  the  real  sensor  sketched  in  equation  2  corresponds  to  a 
projection  of  the  vector  v  +  5v  onto  7£(C-f-6C)  (solid  line  in  figure  2).  Different  manipulation 
matrices  produce  different  approximate  solutions.  For  instance,  the  length  OA  corresponds 
to  the  choice  M  =  (1  0),  i.e.  to  neglecting  the  second  measurement;  OC  corresponds  to 
M  =  (0  1),  neglecting  the  first  one;  OB  represents  the  least-squares  approximation  obtained 
with  M  =  CT .  The  “true”  solution,  of  course  unknowable,  can  now  be  represented  by  the 
length  OY.  The  total  error  Sp  is  the  sum  of  a  projection  error  (e.g.,  OB  —  OY  in  the  case 
of  least  squares),  and  of  an  algorithmic  error,  which  propagates  in  solving  equation  5  (i.e.  in 
finding  the  length  OB). 

The  sensor  design,  i.e.  the  choice  of  the  elements  of  C,  place  the  subspace  V(C)  in 
5?2.  By  varying  this  choice,  and  keeping  6v  and  6C  constant,  the  projection  error  does  not 
change,  since  this  variation  amounts  to  a  rigid  rotation  of  the  graph  of  figure  2  about  the 
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point  0 .  Therefore,  it  seems  reasonable  to  apply  design  criteria  that  tend  to  optimize  the 
accuracy  in  evaluating  p  ,  the  best  affordable  approximation  of  the  “true”  load  with  a  given 
technological  implementation. 

Applying  equation  3  to  equation  5,  we  obtain: 


—  ( £mv  €Tnc  )ifp(MC), 


(7) 


where  er 


||f(Mv) 

!|Mv|j 


[  _  |[*5(MC)|! 

’  tmC  ~  :|MC||  ' 


In  order  to  correctly  compare  sensors  that  employ  different  solution  methods  for  equa¬ 
tion  5,  or  even  sensors  built  with  different  numbers  of  basic  transducers  but  comparable 
technology,  we  need  to  express  equation  7  explicitly  in  terms  of  relative  source  errors  ev, 
tc.  By  “comparable  technology”  it  is  meant  here  that  all  basic  transducers  have  the  same 
statistical  properties  in  terms  of  measurement  and  calibration  accuracy,  such  that  e„  and  ec 
are  not  varied  by  adding  or  removing  transducers  3. 

The  term  emc  can  be  rearranged  as 


||(M  +  <5M)(C  +  8C)  —  MCjj  _ 
IIMCII 

!|M  6C  +SM  C||  „  l|M|[  ||C|| 

||MC||  -  ||MC||  '£m 


where  's  the  relative  error  the  matrix  M  is  known  with.  As  for  the  term  em„,  it 

holds: 


|(M  +  £M)(v  4-  £v)  —  Mv| 

w 

|M  6v  -f  £M  v 


l|Mv[| 


< 


|M| 


iMvl 


(^m  +  €„). 


(8) 


In  order  to  give  a  proper  upper  limit  to  the  last  term  of  equation  8,  which  propagates 
source  errors  into  the  computed  measurements  Mv,  it  is  expedient  to  consider  only  vector 
and  matrix  2-norms  (i.e.  the  Euclidean  and  maximum  singular  value  norm,  respectively). 
It  should  be  pointed  out  that,  although  this  choice  was  not  specified  before  for  the  sake  of 
generality,  the  use  of  2-norms  is  always  advisable  in  problems  like  the  present  one.  In  fact 

3There  may  actually  he  an  improvement  of  eT  by  adding  more  transducers  with  the  same  individual 
accuracy,  in  the  sense  that  more  measurements  of  the  same  quantities  would  be  available.  The  same  effect 
could  be  obtained  on  the  other  hand  by  repeating  the  measurement  of  some  or  all  transducers  and  averaging 
before  solving  the  equations.  We  assume  for  brevity’s  sake  that  both  those  effects  are  accounted  for  in 
evaluating  e„ . 
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these  norms  have  the  property  of  being  invariant  with  respect  to  orthogonal  transformations, 
due  e.g.  to  rigid  1  stations  of  the  reference  frame.  On  the  other  hand,  being  the  dimensions 
of  vectors  and  matrices  involved  in  force  sensing  small,  the  computational  load  is  not  much 
larger  for  2-norms  than  for  other  possible  norms.  In  the  following  only  2-norms  will  be 
considered,  and  explicit  notation  will  be  omitted. 

The  ratio 


^mv 
^  m  “H 


< 


liMIIIMI  ||M|| 
»Mv|| 


(9) 


might  appear  unbound  at  varying  v.  In  particular,  for  v  lying  in  the  null  space  of  M,  j\J  (M). 
infinite  amplification  of  source  errors  would  occur.  This  corresponds  to  a  choice  of  M  that 
neglects  every  meaningful  information  from  the  sensor.  More  precisely,  an  upper  bound  to 
the  ratio  9  can  be  obtained  as  follows.  Consider  the  orthogonal  decomposition  of  the  vector 
v  =  Vp-f  vh,  with  vh  €  S/(M),  and  vp  £  72.(MT).  Using  equation  1,  we  can  write  vh  =  C  ph. 
This  would  not  be  possible  if  v/,  ^  bu,  in  such  case  V/,  vould  not  correspond  to  a. 

physically  meaningful  load  configuration,  so  that  it  must  be  £  ^I(C).  On  the  other  hand, 
if  we  apply  equation  5  we  have  that 


Mvt  =  0  =  MCph. 


This  equation,  along  with  the  previous  assumption  det(MC)  ^  0,  leads  to  the  conclusion 
ph  =  vh  =  0. 

Since  for  what  has  been  shown  the  existence  of  components  of  v  belonging  to  Af(M.)  can 
be  excluded,  equation  9  can  be  modified  to 

<  l|M||  =  IIMJI 

Cm  +  C«  infVp  /‘"in 

where  nmin  indicates  the  smallest  non-zero  singular  value  of  M  (i.e.,  the  square  root  of  the 
smallest  non-zero  eigenvalue  of  MrM).  It  can  be  easily  verified  that 

M  =  ||MMt||  =  iV(M)  =  iVM, 

ft  min 

where  the  obvious  generalization  of  the  definition  of  2-norm-based  condition  number  to 
rectangular  matrices  is  introduced  [21]. 

Finally,  equation  7  can  be  restated  as 


£p  < 


IMil 


IMCI 


-{(■c  +  Cm)  +  Nm((v  +  fm)  >  ATp(MC)  =  F(C,M,ev,ec,e„ 


(10) 
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This  formula  is  the  most  general  result  of  this  paper.  The  design  objective  function  F(-) 
is  defined  as  a  function  of  the  compliance  matrix  C,  of  the  adopted  manipulation  M,  and 
of  the  relative  source  errors  ec,  ev,  and  em.  Note  that  C  and  M  are  in  turn  functions  of  the 
design  variables  x.  Most  often  tm  can  be  expressed  in  terms  of  C  and  ec,  so  that  it  will  be 
dropped  from  the  arguments  of  F(-).  Given  the  values  of  relative  source  errors  ev,  ec  that 
are  expected  with  the  adopted  degree  of  technological  sophistication,  the  accuracy  of  force 
sensors  differing  in  design,  solution  method  and  number  of  transducers  (provided  that  all  the 
transducers  have  similar  statistical  properties)  can  be  compared  by  means  of  equation  10. 

Two  particular  cases  are  very  important  in  practice,  and  are  briefly  discussed  in  the 
following  subsections. 


5.1  Sensors  Solved  With  Least  Squares  Methods 

A  general  least-squares  approximation  of  the  inconsistent  set  of  equations  of  a  redundant 
force  sensor  can  be  obtained  if  the  choice  M  =  WCT  is  made  in  equation  5,  where  W  is 
a  diagonal  matrix  of  weights  corresponding  to  different  accuracies  of  basic  measurements. 
Equation  10  can  be  applied  to  estimate  error  propagation  in  such  a  sensor.  However,  since 
a  correct  weighting  is  rarely  achievable,  the  regular  least-squares  solution  with  M  =  C2  is 
the  only  one  currently  adopted  in  practice  (see  e.g.  [15]). 

Applying  equation  10  and  considering  that,  for  2-norms,  it  holds  ||Cr||  =  |jC]|  and 
||CTC||  =  ||C||2,  we  have: 


eP  <  {[2  +  Nc]  ec  +  Ncev}  KP(CT C),  (11) 

where 

/V2 

*’(cc)  =  rr2fe**-  <12> 

Equation  11  and  equation  12  show  how  the  dependence  of  sensor  accuracy  on  the  condi¬ 
tion  number  is  increased  for  sensors  using  least-squares  approximation  4.  On  the  other  hand, 
redundant  transducers  can  reduce  the  condition  number  Nc-  In  summary,  the  effectiveness 
of  redundant  transducers  should  be  checked  in  the  specific  case,  as  they  may  result  either 

4  It  can  be  noted  that  equation  11  and  equation  12  suggest  a  more  pessimistic  estimate  of  error  propagation 
in  least  square  problems  than  is  currently  held  in  numerical  stability  analysis  [21].  Our  bound  on  error 
propagation  is  in  effect  quite  conservative;  propagation  problems  with  redundant  sensors  can  be  relaxed  using 
more  sophisticated  algorithms  for  the  inversion  of  the  normal  equation  5  with  M  =  CT,  e.g.  Householder 
or  SVD  methods  [23].  More  precise  bounds  would  however  involve  an  a  posteriori  evaluation  of  errors, 
which  is  not  clearly  feasible  in  our  problem.  Relations  11  and  12  provide  a  concise  objective  function  for 
optimization  procedures,  whose  degree  of  conservatism  is  probably  of  the  same  order  for  different  designs 
put  in  comparison. 
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beneficial,  or  useless,  or  even  harmful  for  sensor  precision  (this  point  will  be  reconsidered 
later). 


5.2  Sensors  With  Structured  Matrices 

Some  force  sensor  designs  produce  compliance  matrices  whose  elements  obey  to  „uch  a  pat¬ 
tern  that  the  solution  of  the  associated  linear  system  can  be  performed  very  easily.  Well 
contrived  sensor  designs  in  fact  can  realize  diagonal  or  near  diagonal  MC  matrices,  thus  de¬ 
coupling  the  effect  of  each  load  component  on  the  manipulated  measurement  vector  Mv.  A 
good  example  is  described  in  [5].  Although  this  concept  is  not  limited  to  redundant  sensors, 
it  is  easier  to  design  decoupled  sensors  when  m  >  n. 

Good  decoupling  was  a  strict  requirement  for  force  sensors  designed  before  the  advent 
of  present  cheap  and  powerful  electronic  computers.  In  fact,  all  the  data  manipulations 
needed  for  such  sensors  could  be  carried  out  by  simple  analog  summers  and  multipliers, 
most  often  realized  by  arranging  the  Wheatstone  bridges  and  the  amplifier  gains  used  for 
measurement  in  a  peculiar  combination.  In  some  cases,  for  instance  when  requirements 
on  sensor  bandwidth  dominate  over  accuracy  needs,  such  designs  can  be  still  today  used 
effectively  [20]. 

The  M  manipulation  matrix  for  decoupled  sensors  consists  of  a  simple  pattern  of  integer 
numbers,  so  that  we  can  assume  <fM  ~  0.  Equation  10  can  therefore  be  simplified  to  the 
form: 


«,  <  {  X,(MC).  (13) 

Some  interesting  conclusions  can  be  drawn  from  equation  13.  It  can  be  noted  for  instance 
that  a  rescaling  of  the  sensor  equations,  tending  to  equilibrate  them  by  multiplying  the 
characteristic  matrix  by  a  diagonal  weighting  matrix  M  (which  obtains  a  lower  condition 
number  for  MC),  does  not  produce  any  enhancement  of  sensor  accuracy.  This  fact  has 
obvious  physical  meaning,  and  is  in  effect  equivalent  to  Bauer’s  theorem,  well  known  in 
numerical  analysis  5  [22], 

Furthermore,  if  sensors  are  considered  as  particular  redundant  sensors,  we  expect  that 
the  objective  function  calculated  with  equation  3  should  agree  with  the  one  resulting  from 
equation  10,  in  its  specialized  form  13.  This  can  be  in  fact  easily  verified  by  substituting 
M  =  (lnxn|0„x(m-n))  in  equation  13. 

Finally,  equation  13  allows  us  to  address  the  problem  raised  in  section  4,  i.e.  that  the  two 


matrices  Cj 


1  0 
0  1 


and  C2 


(  1  0 
1  0 

\0  1 


have  different  condition  numbers,  although 


sNote  that  this  scaling  is  purely  algebraic,  as  opposed  to  a  weighting  based  on  statistical  properties  of 
measurements,  which  could  in  fact  increase  the  sensor  accuracy,  affecting  the  source  error  e„ 
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the  same  information  seems,  so  to  say,  to  be  stored  in  them.  The  point,  which  has  general 
validity,  is  that  it  makes  no  sense  to  consider  the  objective  function  value  associated  with  the 
sensor  C2,  without  specifying  how  the  equations  are  going  to  be  solved,  i.e.  which  matrix  M 
is  to  be  used.  So,  if  we  decide  to  use  least  squares  methods,  application  of  equation  11  with 
ee  =  0  gives  the  objective  function  value  F( C2,C^’,e0,ec)  =  2y/2ev,  which  is  much  worse 

than  F(C1}I,  Cv,  ec)  =  e„.  Yet,  a  manipulation  M  =  (  ^  !?  ^  can  be  applied,  such  that 


F(C2,M,ev,ec)  —  ev. 

Among  the  factors  that  tend  to  make  the  decoupling  approach  obsolete  in  force  sensor 
design,  one  is  the  complexity  of  mechanical  structures  necessary  to  implement  it.  Design 
complexity  conflicts  with  requirements  on  size,  cost,  and  reliability  of  the  sensor.  Further¬ 
more,  and  more  subtly,  the  complexity  of  structures  (where  for  instance  stress  concentrations 
are  more  or  less  deliberately  introduced),  makes  the  linear  model  on  which  the  sensor  is  based 
even  more  arbitrary.  Besides  that,  it  should  be  pointed  out  that  force  sensors  have  structured 
compliance  matrices  only  as  far  as  they  are  on  the  designer  drawing  table.  No  calibration 
technique  is  allowed  by  decoupling  design  to  correct  implementation  dependent  errors  on  C, 
so  that  larger  errors  ec  must  be  expected  (typically  an  order  of  magnitude  larger  than  in 
calibrated  sensors,  in  the  author’s  experience). 


6  Accuracy  And  Adding  More  Transducers 


Equation  10  identifies  an  objective  function  capable  of  leading  an  automatic  design  proce¬ 
dure  in  the  choice  of  optimal  solutions.  Nevertheless,  a  better  understanding  of  transducer 
redundancy  can  be  gained  by  analyzing  how  Nc  is  varied  by  adding  new  transducers. 

If  C  €  3?"*xn  is  the  compliance  matrix  of  a  given  (stage  in  the  design  of  a)  sensor, 
an  additional  set  of  d  basic  transducers  (henceforth  called  gauges)  will  produce  the  new 

characteristic  matrix  A  =  ^  ^  ^ ,  where  D  6  !R*xn,  A  €  3^m+<^xn.  The  condition  number 

of  the  new  linear  system,  according  to  the  choice  of  2-norms,  is  the  square  root  of  the  ratio 
between  the  largest  and  the  smelliest  eigenvalues  of  ATA.  Considering  ArA  as  the  sum  of 
two  symmetric  matrices, 

ArA  =  (  CT  |  DT  )  =  CrC  +  DtD, 


and  applying  the  minimax  theorem  [23 j,  we  have: 

max^i,n{7,2  +  2  7i  +  % 

mini=1,n{7/  +  £2_i+1}  ~  A  ~  72  +  S2n  ’ 


(14) 


14 


where  (7i,  72,  ■  •  • ,  7n)  and  (tfj,  S2, . . . ,  6n)  are  the  singular  values  of  C  and  D  in  non  increasing 
order,  respectively.  This  relationship  gives  bounds  for  condition  number  variations  as  a 
consequence  of  the  addition  of  new  gauges,  but  it  is  insufficient  to  determine  whether  the 
variations  are  favorable  or  not. 

More  detailed  information  can  be  obtained  if  extra  gauges  are  added  one  at  a  time.  The 
matrix  D  in  this  case  is  reduced  to  a  row-matrix  dr  €  3ilxn,  and  DrD  —  drd  G  3?nxn  has 
singular  values  Si  =  (drd)^2,  S2  =  63  =  •  •  •  =  6n  =  0.  Any  set  of  n  orthogonal  vectors  which 
includes  d  forms  an  eigenvector  set  for  DrD.  Equation  14  is  modified  in  this  case  as 

max{72,72  +  S2}  „2  7i  + 

min{7n-n  7^  +  $1}  “  A  ~  ll  +  Q  ' 


If  d  is  chosen  parallel  to  the  eigenvector  of  CTC  corresponding  to  the  smallest  eigenvalue 
72,  CrC  and  DTD  share  the  same  eigenvector  set.  Therefore,  also  ATA  =  CTC  -f  DrD 
shares  the  same  eigenvectors,  and  the  eigenvalues  of  ArA  are  equal  to  those  of  CrC  except 
for  one,  a\  =  72  +  dTd,  with  1  <  s  <  n. 

Corresponding  to  this  choice  of  the  direction  of  d,  we  have  the  new  condition  number: 


_  max^1|n  {7t2,72  +  drd} 

"  +  drd}  “ 


(15) 


According  to  this  result,  a  new  transducer  should  be  placed  on  the  sensor  in  a  position 
such  that  its  output  um+1  is  related  to  the  load  p  as  um+i  =  dTp,  with  d  parallel  to  the 
singular  vector  corresponding  to  the  minimum  singular  value  of  the  old  characteristic  matrix, 
and  such  that  drd  =  7J  —  72.  By  satisfying  perfectly  these  conditions  with  n  —  1  additional 
transducers,  it  is  theoretically  possible  to  modify  any  given  design  so  as  to  obtain  unitary 
condition  number.  In  practice,  this  goal  may  be  impossible  to  achieve,  due  to  the  fact  that  an 
arbitrary  output  law  can  not  always  be  obtained  from  a  transducer.  However,  these  design 
suggestions,  along  with  equation  10,  can  be  useful  to  decide  whether,  and  how,  to  correct 
sensor  designs  by  adding  transducers. 


7  Applications 

In  order  to  illustrate  some  aspects  of  the  methods  above  outlined,  two  applications  of  op¬ 
timal  force  sensor  design  are  presented  in  this  section.  The  first  one  is  extremely  simple, 
and  intended  only  to  illustrative  purposes.  The  second  example  reflects  a  real  application 
problem,  regarding  the  design  of  miniaturized  force  sensors  designed  to  fit  robotic  fingertips. 
The  sensor  is  currently  being  produced  and  successfully  employed  in  several  applications. 
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1 


2 


Figure  3:  A  simple  cantilever  beam  sensor  with  rectangular  cross  section. 

7.1  A  Simple  Two  Axis  Force  Sensor 

An  extremely  simple  two  axis  force  sensor  can  be  realized  as  depicted  in  figure  3.  By  means 
of  strain  measurements  on  the  cantilever  beam  with  rectangular  cross  section  b  x  h  and 
length  Z,  the  unknown  load  components  pa  and  p2  are  sensed.  The  five  components  of  the 
vector  x  of  design  variables  are  defined  as  Xi  =  b,  z2  =  h ,  and  i3,i4,i5  as  the  distances 
of  the  gauges  from  the  tip  of  the  beam.  Three  unilateral  constraints  express  limitations  on 
acceptable  stress  and  displacement  under  the  maximum  nominal  load  pmox,  as  well  as  on 
the  strain  under  the  minimum  nominal  load  pmjn: 

».(*)  =  ^  +  P-±jp-  -  Vadm  <  0, 

-  /-«  <  o,  (16) 

/  \  Pl.min  P2,minx2  .  e  ~  n 

g3(x)  = - —  +  - — -  +Zmtn  <  0, 

where  A  =  iix2  is  the  cross  section  area,  J  =  ijXj/12  its  momentum  of  inertia,  E  is  the 
Young’s  modulus  of  elasticity  of  the  material,  and  <Tadm  its  maximum  admissible  strain: 
fmax  is  the  maximum  acceptable  displacement  under  load,  and  £mm  the  minimum  sensitivity 
required  to  the  sensor.  The  feasibility  region  is  defined  as 
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0  <  Xi  <  /,  i  =  3, 5, 


while  obvious  bounds  are  set  on  xx  =  b  and  x2  =  h  by  constructive  considerations.  Three 
cases  will  be  considered:  a  minimal  sensor  solved  with  full  inversion  or  with  special  solution, 
and  a  redundant  sensor  with  three  basic  transducers. 

7.1.1  Minimal  realization  with  full  inversion 

For  a  minimal  sensor  using  only  two  gauges,  sx  and  S2>  the  strains  are: 


Pi, max 

Pi 

,  ?2, max  ^2 

EA 

Pi, max 

+  EJ 

x3 

Pi, max 

Pi 

P2,maxx2 

EA 

Pi, max 

EJ 

X4 

The  sensor  model  is  v  =  Cp,  where  p  =  (— 21 - 22 — )  is  the  normalized  load  vector, 

r  r  \  Pi  .max  Pi, max  /  ’ 

and  the  compliance  matrix  (modulo  a  constant)  is 


Pi  1ma£  P2,maxx2  xl 
EA  EJ 

Pi  ,mq»  _ P7,max^2xA 

EA  EJ 


In  this  simple  case,  the  optimum  design  can  be  easily  found.  Imposing  Nq 
to  an  orthogonality  condition  for  the  columns  of  C,  cx  and  c2: 


=  I  amounts 


cfc2  =  0, 

lldll  =  INI 

whence  the  design 

_  *^2  Pi, max 

13  “  14  "  lWf  rZZ' 

Given  such  dimensions  the  resulting  error  is  ep  ~  ev  +  ec.  Constrained  minimization 
techniques  should  be  used  if  these  dimensions  do  not  comply  with  constraints  16. 
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7.1.2  Minimal  realization  exploiting  the  structure  of  C 


If  the  peculiar  structure  that  C  has  for  x3  =  £4  =  x  is  exploited  by  using  a  manipulation 
matrix 


M  = 


the  original  equation  is  modified  to 


Mv  = 


v2  +  v2 
Vi  -  v2 


=  MCp  = 


)  Pl  ,r 


EA 

0 


v  P2 .  1 


.IJl 


P, 


EJ 


which  can  be  solved  with  only  n  =  2  flops,  instead  of  the  n2  =  4  otherwise  necessary. 
This  advantage  is  obviously  of  some,  if  little,  importance  only  for  larger  n  than  used  in  this 
example. 

The  following  relationships  can  be  easily  verified: 


|1VL||  =  ||M  x\\  =  ^2;  Nm  =  1; 
I C II  =  V2  max 


I  —  ,max  P2,maxx2x  1 

V  EA  ’  EJ  i  ’ 


||MC||  =  2  max 
iV(MC)  =  max 


\maxx  2*  1 


f  Pl,max  ?2 ,n 

\~EA~  ’  £/  /; 

f  J Pl,max  Ap2rTnaxX2X 

l  ^PZ.ma 


According  to  equation  13,  if  x  = 


JPl.r 


,maxx2x  JPl  .maj 

is  chosen,  we  have 


€p  ^  | "  ||mc||  “ £c  +  ATp(MC)  «  £v  +  ec, 

which  shows  that  this  manipulation  does  not  worsen  error  propagation.  However,  as  al¬ 
ready  mentioned,  in  practical  cases  the  cancellation  of  the  off-diagonal  entries  of  MC  is  not 
expected  to  be  perfect,  so  that  larger  ec  must  be  taken  into  account. 


7.1.3  Redundant  realization 


If  the  three  gauges  s1?  s2,s3  are  placed  on  the  beam,  the  redundant  sensor  compliance  matrix 
results 


Pl  ,ma x  P2  ,mat  x2 

EA  EJ 

Pl  ,ma»  _ P2,maxx7  x4 

EA  EJ 

Pl  .mar  P2,ma»  J 2 x£ 

EA  EJ 


\ 
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By  imposing  Nc  =  1,  we  obtain  the  optimal  design: 


x3  -  +  xs 

X3  +  X4  +  X\ 

and,  applying  equation  11,  we  have 


0 

_3_  Pjm^l 
144  plmax 


£p  £  3ec  +  ev. 

This  result  shows  that  the  introduction  of  the  third  gauge  makes  error  propagation  worse. 
On  the  other  hand,  the  optimality  conditions  above  allow  a  wider  design  flexibility,  which 
may  ue  necessary,  depending  upon  the  constraints  applied  to  the  variable  ranges. 


7.2  A  Six  Axis  Miniaturized  Force  Sensor 

As  mentioned  in  the  introduction,  a  particularly  promising  application  of  force  sensors  is 
intrinsic  contact  sensing  [8].  In  its  most  common  embodiment,  this  method  adopts  six-axis 
force  sensors  integrated  as  close  as  possible  to  the  parts  of  the  arm  whose  surfaces  contact 
the  environment.  This  application  has  critical  size  requirements:  for  instance,  force  sensors 
have  been  designed  to  fit  the  interior  of  fingertips  of  articulated  dextrous  hands  as  small  as  a 
cylinder  18mm  diameter  and  25mm  long.  Furthermore,  contact  sensing  algorithms  use  force 
sensors  outputs  as  input  data,  thus  making  accuracy  requirements  also  very  important. 

The  approach  to  the  design  of  sensors  for  intrinsic  contact  sensing  that  will  be  illustrated 
here  basically  follows  the  guidelines  illustrated  below,  which  we  consider  in  general  advisable 
for  high-performance  sensor  design: 

•  Use  very  simple  mechanical  structures,  whose  behavior  is  as  close  to  linear  as  possible; 

•  Utilize  the  least  necessary  number  of  transducers,  to  avoid  unnecessary  propagation  of 
errors; 

•  Give  the  design  a  convenient  parameterization; 

«  Use  numerical  techniques  to  find  the  best  combination  of  design  variables  which  com¬ 
plies  with  design  constraints. 

Following  this  approach,  a  miniaturized  force  sensor  has  been  designed  [24]  and  built 
in  several  copies,  used  to  sensorize  devices  for  dextrous  manipulation  such  as  the  Salisbury 
hand  [25],  the  “tactile  explorator”  finger  of  the  Centro  “E.Piaggio”  [11],  and  the  whole-hand 
manipulation  system  UB-Hand  II  [3].  The  mechanical  structure  of  this  sensor  is  sketched 
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Figure  4:  A  miniaturized  force/torque  sensor  with  cylindrical  beam  structure. 

in  figure  4:  it  consists  of  an  hollow,  thin-walled  cylinder.  Strain-gauges  are  applied  on  the 
external  surface  of  the  cylinder,  in  a  number  of  six  6.  The  cylinder  dimensions,  the  position 
of  the  gauges  on  its  structure  and  their  orientation  have  been  considered  as  the  variables  to 
be  optimized. 

Figure  4  shows  the  {0,£i,  £2,(^3}  reference  frame  in  which  the  components  of  the  load 
applied  to  the  sensor  extremities  are  described.  The  axis  is  placed  along  the  cylinder  axis. 
The  position  of  the  itf%  gauge  is  uniquely  determined  by  the  cylindrical  coordinates  of  its 
center  point,  C.,  and  0,,  and  by  the  angle  fa  formed  by  the  gauge  axis  with  the  cylinder  axis. 
The  design  variables  are  three  for  each  gauge,  plus  the  cylinder  radius  and  wall  thickness, 
i.e.  20. 

The  simplicity  of  the  mechanical  structure  allows  the  evaluation  of  the  entries  of  the 
compliance  matrix  C  by  means  of  simple  relations  of  elastic  beam  theory  [26].  Note  that  the 
computation  of  C  has  to  be  done  at  each  iteration  step  of  the  numerical  optimizing  routine, 
so  that  avoiding  demanding  techniques  as  the  finite  element  method  is  extremely  expedient. 
Accordingly,  the  strain  corresponding  to  the  ith  gauge  is 

vi  =  ^2  Cx.iVji 
.7=1,6 

9This  number  is  actually  increased  to  7  for  temperature  compensation;  this  point  is  addressed  in  the 
appendix. 
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where  p:  is  the  value  of  the  jth  load  component  normalized  with  respect  to  its  nominal 
maximum  value,  pjiTnax.  Elements  CX<J  are  as  follows: 


Ci,i 

—  Pi, max 

kPn(cos2  4>i  —  i/sin2 

fa); 

Ci,  2 

—  P2,max 

[W/C* cos  ^,-(cos2  <f>i 

—  u  sin2  <t>i)  + 

Wa  sin  &i 

sin  2&]  ; 

Ct>3 

~  p3,max 

[vP/Cj  sin  ^(cos2  <f>i 

—  v  sin2  4>i)  + 

W3  cos  9, 

sin  2 ; 

C;,4 

P  Aim  ax 

Wt  sin  2  fa 

Cit6 

—  Ps,max 

W ;  cos  ^,(cos2  <t>i  — 

v  sin2 

Ci,  6 

~  P6,max 

Wf  sin  #i(cos2  4>i  — 

is  sin2  <pi); 

where  v  is  the  Poisson’s  ratio  for  the  structure  material,  and  the  moduli  W  are  defined  as 

1 

2irRsE' 

2  Wn 

R  ' 

2(1  +  u)Wn; 

(1  +  v)Wn 

R 

where  E  is  Young’s  modulus,  R  is  the  cylinder  radius  and  s  its  wall  thickness.  The  relations 
above  are  valid  in  the  assumption  s  R:  in  this  case,  l/s  is  a  common  factor  of  C%,:,  so 
that  Nc  is  not  affected  by  s.  In  other  words,  the  relative  accuracy  of  the  sensor  does  not 
depend  on  the  wall  thickness;  s  can  be  chosen  independently  of  other  design  variables,  and 
made  such  that  minimum  sensitivity  and  maximum  strain  level  requirements  are  met. 
Other  design  constraints  for  a  sensor  to  fit  a  robotic  hand  fingertip  are  of  the  form 


Wn  = 

Wf  = 
W.  = 
wt  = 


0  >  ICi  -  01  -  Imam, 

Rmax  >  R>2S, 

where  lmax  =  15 mm,  Rmax  =  6 mm  are  typical  values. 

According  to  the  previous  discussion,  an  objective  function  F(R,£i,9i,<j>i)  =  Nc  has  been 
assumed  to  optimize  this  minimal  sensor.  Note  that  for  minimal  sensors,  we  do  not  need  to 
use  explicit  values  of  the  source  errors  and  ec.  The  condition  number  of  the  compliance 
matrix  is  a  complex  function  of  the  19  design  variables  whose  minima  could  not  be  found 
analytically.  A  numerical  algorithm  (a  customized  version  of  Powell’s  and  Broyden’s  methods 
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[22])  has  been  employed.  The  values  of  design  variables  corresponding  to  the  lowest  of  the 
relative  minima  found  in  repeated  trials  (resulting  in  Nc,min  —  3.07),  are  given  in  table  1. 
A  singular  value  decomposition  of  the  corresponding  compliance  matrix  is  reported  in  table 
2,  showing  that  a  3.07  ratio  occurs  between  the  maximum  and  minimum  relative  sensitivity, 
the  latter  corresponding  approximately  to  the  pi  component  of  the  load,  acting  along  the 
cylinder  axis. 


Variable 
Gage  no. 

p(rad) 

0(rad) 

Cl  nun) 

R(mm) 

s(rrun) 

1 

0.3 

-0.7 

1.2 

2 

-0.3 

1.6 

1.7 

3 

0.2 

3.2 

-1.9 

5.4 

0.2 

4 

-0.6 

-1.8 

-1.4 

5 

0.4 

0.3 

-5.9 

6 

0.3 

0.9 

-6.0 

_ 

Table  1:  Optimal  design  variables  for  the  six-axis  miniaturized  force  sensor,  corresponding 
to  the  nominal  load:  px  =  p2  =  p3  =  107V;  p4  =  SON  mm;  p5  =  p6  =  85  Nmm. 


Singular  Structure 

Singular 

7i 

72 

73 

74 

75 

76 

values 

138.4 

129.6 

89.4 

64.3 

55.9 

45.0 

Singular 

«1 

u2 

u3 

U4 

«5 

U6 

vectors 

0.07 

0.08 

0.00 

-0.05 

-0.27 

0.95 

-0.04 

0.32 

0.83 

-0.03 

-.17 

-0.04 

0.05 

-0.36 

0.10 

0.81 

-0.44 

-0.05 

0.08 

-0.18 

0.31 

0.34 

0.82 

0.26 

0.02 

0.83 

-0.28 

0.46 

0.10 

-0.02 

0.91 

0.16 

0.34 

-0.10 

-0.11 

-0.12 

Table  2:  Singular  values  and  corresponding  singular  vectors  of  the  sensor  compliance  matrix 
for  the  optimal  set  of  design  variables  listed  in  table  1. 
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Figure  5:  Sub-optimal  arrangement  of  gauges  on  the  miniaturized  cylindrical  sensor. 


7.2.1  Sub-optimal  design 

An  interesting  sub-optimal  configuration  of  the  same  force  sensor  is  depicted  in  figure  5.  The 
centers  of  the  gauges  are  placed  on  two  circles,  on  the  plane  (j  =  0  for  gauges  1  through 
4.  and  on  the  plane  £i  =  —d  for  gauges  5  and  6.  Angles  are  6 ,  =  i~ ,  and  <f>,  =  (  — 

1  £  i  <  6.  Design  variables  are  thus  reduced  to  three,  R,  d ,  and<z>.  The  compliance  matrix 
given  by  such  design  has  the  following  structure: 

/  a  0  -b  e  0  /  \ 

a  b  0  -e  /  0  } 

_  _  a  0  6  e  0  — / 

a  -b  0  -e  — /  0 

a  c  -b  e  0  / 

\  a  b  —  c  —  e  /  Oy 

An  almost-decoupled  sensor  can  then  be  obtained  even  with  this  extremely  simple  design, 
by  using  a  manipulation  matrix 
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M  = 


/  1 
-1 
0 
1 
0 

V  i 


i 

o 

i 

-l 

i 

0 


1  1  0  0  \ 

0  0  10 

0  0  0  -1 

1-10  0 
0-100 
-1  0  0  0  / 


An  extremely  fast 


solution  of  the  sensor  equations  is  possible  with 


this 


arrangement: 


1  VX 

+  v2 

+ 

v-i 

t  tq  ' 

1  4g 

0 

0 

0 

0 

0  ^ 

Vs 

-  Vi 

0 

c 

0 

0 

0 

0 

v6 

—  v2 

0 

0 

c 

0 

0 

0 

t>i 

-  v2 

+ 

V3 

-  v4 

— 

0 

0 

0 

4e 

0 

0 

P 

v2 

-  VA 

0 

2b 

0 

0 

2/ 

0 

\  vi 

-  V3 

) 

l  0 

0 

-2b 

0 

0 

2/  / 

The  optimal  combination  of  design  variables  with  sub-optimality  constraints  is  reported 
in  table  3.  The  corresponding  objective  function  minimum  is  at  Nc  =  6.6.  As  expected,  the 
sub-optimal  sensor  is  less  accurate  than  the  optimal  one.  The  advantage  of  this  design  is 
that  it  is  possible  to  choose  either  a  precise  solution  algorithm  (full  inversion  of  C)  or  a  fast 
one.  Perhaps  more  importantly  from  a  practical  point  of  view,  the  symmetries  of  the  design 
render  its  fabrication  much  easier  than  the  fully  optimized  one.  It  is  finally  noted  that  the 
objective  function  is  fairly  smooth  in  the  neighborhood  of  its  minimum,  which  fact  renders 
the  design  quite  robust  to  fabrication  inaccuracies. 


Variable 
Gage  no. 

0(rad) 

C(mm) 

i2(mm) 

s(mm) 

1 

0.3 

0 

-7.9 

2 

-0.3 

7t/2 

-7.9 

3 

0.3 

7T 

-7.9 

5.6 

0.2 

4 

-0.3 

3x/2 

-7.9 

5 

0.3 

0 

4.1 

6 

-0.3 

t/2 

4.1 

Table  3:  Sub-optimal  design  variables  for  the  six-axis  miniaturized  force  sensor,  correspond¬ 
ing  to  the  nominal  load:  P\  =  p2  —  Pi  =  10iV;  p4  =  60Nmm ;  ps  =  p6  =  85 Nmm. 
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Table  4:  Singular  values  and  corresponding  singular  vectors  of  the  sensor  compliance  matrix 
for  the  sub-optimal  set  of  design  variables  listed  in  table  3. 


7.2.2  Some  Experimental  Data 

As  it  has  been  pointed  out  several  times  already,  the  ultimate  performance  of  a  force  sensor 
depends  on  design  quality  as  well  as  on  technological  factors.  It  is  not  therefore  easy  to  assess 
design  quality  by  experimental  means,  as  an  extremely  large  statistical  basis,  and  a  tight 
control  of  experimental  conditions  (including  details  such  as  quality  of  gauges,  of  bonding 
agents,  of  electronic  components  and  so  on),  would  be  required.  Although  such  exhaustive 
verification  can  not  be  provided  here,  some  numerical  values  obtained  from  miniaturized 
multi-axis  sensors  will  be  given  for  generic  reference. 

An  aluminum  (2024  T4)  sensor  designed  according  to  the  sub-optimal  scheme,  with  wall 
thickness  0.8mm  and  using  foil-gauge  components,  has  been  connected  to  instrumentation 
amplifiers  and  to  an  A/D  port  of  a  personal  computer.  The  maximum  nominal  load  was 
30N  for  force  components  and  150Nmm  for  moment  components.  An  estimate  of  the  relative 
source  errors  is  ty  «  ec  ~  1%.  As  a  figure  of  the  global  accuracy  of  a  sensor,  the  cross-talk 
matrix  X  is  sometimes  used,  whose  Xij  entry  is  defined  as  the  ratio  between  the  measured 
/alue  of  the  ith  component  of  the  nominal  load,  Pi,  and  the  actual  value  pj  (for  an  ideal 
sensor,  X  =  I).  The  cross-talk  matrix  obtained  for  the  sensor  (employing  a  full  inversion 
solution  scheme)  is  as  follows: 
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/  1.02 

-0.01 

0.001 

-0.01 

0.01 

0.01 

0.01 

1.01 

-0.03 

-0.01 

0.00 

0.00 

0.01 

0.01 

0.99 

0.05 

0.00 

0.01 

0.01 

-0.02 

0.00 

0.98 

-0.01 

-0.02 

0.00 

-0.01 

-0.02 

-0.03 

1.02 

0.01 

V  -o.oi 

-0.01 

-0.01 

0.04 

0.04 

0.99 

) 

Assuming  the  maximum  norm  of  the  columns  of  X  as  a  measure  of  the  relative  error, 
tp  =  7.5%  can  be  obtained  by  this  matrix,  which  is  lower  than  the  expected  bound  ep  < 
Nc  ~  13%.  This  result  is  believed  to  depend  partly  on  the  conservative  nature  of  the 
condition  number  criterion,  partly  on  the  fact  that  the  “worst  possible  combination”  of  load 
components  is  not  reflected  in  the  cross-talk  matrix. 


8  Conclusions 

It  has  been  shown  that  the  design  of  a  multi-axis  force  sensor  can  be  optimized  with  respect 
to  its  accuracy  by  minimizing  an  objective  function  of  the  design  variables.  The  proposed 
function  of  merit  (equation  10)  is  general  enough  to  allow  the  choice  of  the  optimal  design  in 
a  broad  range  of  possible  solutions.  In  particular,  the  choice  of  the  optimal  number  of  basic 
transducers  to  be  used  in  a  multi-component  sensor  is  addressed.  For  sensors  employing  as 
many  basic  transducers  as  the  unknown  load  components  (minimal  sensors),  the  objective 
function  reduces  to  the  condition  number  of  the  sensor  compliance  matrix  (equation  3). 

An  often  asked  question  about  multicomponent  sensing  is:  “How  many  basic  transducers 
should  be  utilized?  Are  redundant  transducers  useful,  or  useless,  or  are  they  harmful  to  the 
overal  sensor  performance?”  A  clear-cut  answer  to  this  question  is  not  possible,  since  in  the 
general  case  the  direct  evaluation  of  equation  10  is  necessary  to  assess  the  results.  However, 
a  general  guideline  that  results  from  the  discussion  of  equation  10  is  that,  whenever  the 
optimal  configuration  of  a  minimal  sensor  complies  with  the  constraints,  adding  redundant 
sensors  is  unlikely  to  improve  sensor  accuracy. 

Finally,  we  would  like  to  note  that,  although  multi- axis  force  sensors  have  been  explic¬ 
itly  considered  in  this  paper,  the  method  proposed  is  more  general,  and  can  be  directly 
extended  to  any  sensor  of  multiple  physical  quantities  based  on  a  number  of  linearly  related 
measurements. 
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Appendix:  Compensation  of  Thermal  Effects 

The  various  effects  of  temperature  variations  are  among  the  most  important  error  sources 
for  basic  strain  transducers.  For  resistive  strain-gauges,  such  effects  can  be  grouped,  in 
decreasing  order  of  importance,  as  follows: 

•  variation  of  the  intrinsic  (non  stressed)  gauge  resistance; 

•  thermal  expansion  of  the  sensor  structure; 

•  gauge-factor  variation; 

•  variation  of  the  elasticity  modulus  of  the  structure  material. 

While  the  first  two  phenomena  cause  errors  in  measurements  vz,  the  latter  two  modify 
the  compliance  matrix  elements,  entering  the  measurement  output  in  a  nonlinear  manner. 
Strain  gauges  exhibiting  a  curve  of  intrinsic  resistance  vs.  temperature  that  approximately 
matches  (at  least  in  some  temperature  range)  the  thermal  expansion  vs.  temperature  char¬ 
acteristic  curve  of  the  material,  are  currently  available.  Such  gauges  perform  a  first,  rough 
compensation  of  thermal  effects. 

Possible  approaches  to  finer  compensation  are  listed  as  follows,  in  increasing  order  of 
accuracy: 

•  To  regard  the  effects  of  temperature  variations  AT  as  a  component  of  the  source  errors. 
The  counteraction  is  to  minimize  error  propagation; 

•  To  regard  nonlinear  effects  of  temperature  as  source  errors,  an  assume  that  the  remain¬ 
ing  effects  are  functions  of  the  temperature  only,  independent  of  the  individual  gauge 
characteristics,  position,  orientation  etc.  Since  the  effects  of  temperature  are  common 
to  every  gauge  on  the  sensor,  a  dummy  gauge  applied  on  a  stiff  part  of  the  structure, 
and  subject  to  the  same  temperature  variations,  can  be  taken  as  a  reference. 

•  To  consider  nonlinear  effects  as  source  errors,  and  model  the  remaining  effects  of  tem¬ 
perature  T  with  suitable  individual  functions  r,(T).  These  functions  are  typically 
approximated  by  polynomials  of  degree  three  or  higher.  Starting  from  data  normally 
provided  by  the  gauge  supplier,  the  width  wt,(T)  of  a  band  within  which  the  r ;  of  every 
gauge  in  the  batch  employed  lies,  can  be  obtained  (see  figure  6. a).  If  two  curves  ri(T) 
and  ru(T)  are  considered  such  that,  for  any  T  and  any  i,  r;(T)  <  r,(T)  <  ru(T),  we 
have: 
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Figure  6:  Typical  shape  of  foil  gauge  resistance  vs.  temperature  curves:  a)  absolute;  b) 
relative  to  an  offset  measurement. 


ru(T)  -  ri(T)  <  wbru(T). 

If  moreover  an  offset  measurement  is  taken  at  T  —  T0,  the  curves  Ar,(T  —  T0)  = 
At\(AT),  can  be  represented  as  in  figure  6.b.  Let  Arm(AT)  =  Ar‘+-r-*  be  the  median 
curve  of  the  band,  so  that  any  curve  A 7%  can  be  approximated  by  A ^  Ct,,Arm(ST), 
where  cti,  is  a  constant  characteristic  of  the  ith  gauge.  An  upper  bound  for  the  approx¬ 
imation  error  is  wb.  The  effects  of  the  temperature  on  a  single  measurement  can  be 
expressed  as 


Vi  =  c,p  +  ct,,Arm(  AT), 

where  C,  is  the  ith  row  of  C.  Equation  10  can  therefore  be  rewritten  as 


v  = 


=  C'p'. 
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The  elements  cti,  of  the  column  vector  ct  can  be  calibrated  in  controlled  temperature 
conditions.  The  compensated  sensor  requires  the  solution  of  a  linear  system  of  m 
equations  in  7  unknowns.  The  seventh  unknown,  Arm,  is  actually  of  little  practical 
interest,  and  its  calculation  can  be  superseded,  as  long  as  the  whole  calibration  matrix 
has  been  used  to  evaluate  the  non-thermal  components  of  load.  The  optimization 
procedures  discussed  in  this  paper  are  still  valid  for  sensor  adopting  this  compensation 
technique,  if  applied  to  the  new  characteristic  matrix  C'. 
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